Reply to Lindgren and Salomonson 



R. K. Nesbet 

IBM Almaden Research Center, 650 Harry Road, San Jose, CA 95120-6099, USA 

(February 2, 2008) 



For Phys.Rev.A 

In the accompanying Comment [Phys. Rev. A 67, 056501 
(2003)], I. Lindgren and S. Salomonson claim to prove for the 
Kohn-Sham kinetic energy functional of ground state electron 
density that a Frechet functional derivative exists, equivalent 
to a multiplicative local potential function. If true, this re- 
sult would imply an exact Thomas-Fermi theory for ground 
states of noninteracting electrons. However, such a theory 
is not consistent with the exclusion principle for more than 
one electron of each spin. The simplest counterexample is the 
lowest triplet state of a noninteracting two-electron atom. If 
only the total electron density were normalized, as in Thomas- 
Fermi theory, the lowest state would collapse into a doubly- 
occupied Is spin-orbital. Two independent parameters ei s 
and £2s are required to maintain independent subshell nor- 
malization. The argument presented by these authors is dis- 
cussed in the light of this unphysical implication. 

The underlying issue here is whether Hohenberg-Kohn 
theory [1] implies an exact ab initio Thomas-Fermi the- 
ory (TFT) for ground states. For noninteracting elec- 
trons, the question reduces to the existence of a Frechet 
functional derivative [2] of the Kohn-Sham kinetic energy 
functional [3] , equivalent to a multiplicative local poten- 
tial function. If such a functional derivative did exist, it 
could be expressed as 

X"=M-«(r), (1) 
Sp 

which simply states the TFT Euler-Lagrange equation. 
Constant p is determined by the normalization condition 
J p = N for N electrons. In the context of density func- 
tional theory (DFT), p = J2i n iPi = J2i n i4>i4>ii a sum of 
orbital subshell densities. Density variations generated 
by unconstrained independent variations in the orbital 
Hilbert space are necessary and sufficient to determine 
the Euler-Lagrange equations. Spin indices and sums are 
assumed but suppressed in the notation here. For nonin- 
teracting electrons, the correct Euler-Lagrange equations 
for Fermi-Dirac electrons, as used in Kohn-Sham theory, 
are a set of Schrodinger equations, 

t<j> i (r) = {e i -v(r)}<p i (r), (2) 

where t = — ^V 2 , for orthogonal but independently nor- 
malized orbital functions 4>i(r). Independent orbital nor- 
malization requires independent parameters e,, and or- 
bital orthogonality implies the exclusion principle. The 
restriction to one free parameter in Eq.(l) means that the 
exclusion principle cannot be enforced unless all energy 



eigenvalues u are equal. Hence, in a compact system, 
these equations cannot be reconciled for more than one 
electron of each spin. The simplest specific example is 
the lowest ls2s 3 S state of a two-electron system. TFT 
implies bosonic condensation with two electrons in the Is 
spin-orbital. The argument of Lindgren and Salomonson 
[4] (LS) must fail, because it implies that TFT is an exact 
theory for noninteracting electrons. The proposed argu- 
ment also disagrees with DFT, in which the Kohn-Sham 
equations for noninteracting electrons are just Eqs.(2), 
valid for atomic shell structure and consistent with the 
exclusion principle. The TFT and DFT Euler-Lagrange 
equations are inconsistent in the ls2s 3 S example, be- 
cause two independent constants cannot be deduced from 
one. This conclusion has been disputed and defended in 
recent publications [5-7]. 

LS put forward two variational derivations that set 
out to determine the density functional derivative of the 
Kohn-Sham kinetic energy functional T s [p], constructed 
as the ground-state limit of the Schrodinger orbital func- 
tional T[{^>t}] = J2i n>i(i\t\i). In Section IIA, they con- 
sider an N-electron wave function For noninteracting 
electrons, the Schrodinger equation is separable, and the 
wave function can be expressed as a serial product of or- 
bital functions, W = IL</>i(ri). Fermi-Dirac statistics are 
not implied, but must be imposed by explicit orthogonal- 
ization and by independent normalization of the partial 
densities pi = 4>*4>i- In terms of partial densities, LS 
Eq.(18) should be ST — JJ2t n i( e i ~ v )&Pi- Instead, LS 
extend the functional derivative of Eq.(22) to unnormal- 
ized densities by making the implicit assumption that 
the undetermined constant has no orbital subshell struc- 
ture, which is inconsistent with the orbital Schrodinger 
equations. This point can be examined in detail by con- 
sidering specific wave functions, as in LS Section IIB. 

In Section IIB, LS propose a variant of standard 
Hartree-Fock theory in which the product of orbital over- 
lap integrals is retained, in contrast to the usual sim- 
plified argument based on the orbital functional T. If 
done correctly, this variant leads to standard Hartree- 
Fock equations, equivalent for noninteracting electrons 
to Kohn-Sham equations, in which the kinetic energy 
operator t takes its usual form. The problem here is 
that all terms in the total energy expression for a Slater 
determinant must be considered together. To enforce or- 
bital orthonormality, as required by the definition of a 
Slater determinant, a matrix of Lagrange multipliers 
is required. It is customary (for closed-shell systems) to 
assume a diagonalized canonical representation, leaving 
the diagonal elements ti to be determined as eigenval- 
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ues of the canonical Hartree-Fock equations. These are 
Lagrange multipliers for independent normalization con- 
straints = 1. Derivations using the simplified or- 
bital functional T are standard in applications of DFT 
and are consistent with Hartree-Fock theory. Products of 
the overlap factors are unity for variational solutions and 
can consistently be omitted. If retained in the deriva- 
tion, overlap variations are multiplied by a factor that 
vanishes. An appropriate variational functional for non- 
interacting electrons is 



^2m{i\i+v{r) - ei\i)ILjjti{j\j). 



(3) 



This is derived directly for a serial product wave func- 
tion, with Lagrange multipliers for independent orbital 
normalization. Orbital orthogonality is implied by the 
Euler-Lagrange equations or by symmetry. The first fac- 
tor in each term here vanishes in stationary states, so 
all variations of the overlap product drop out, and varia- 
tions of the first factor are multiplied by unity. Written 
out in full, the correct differential for variations about a 
stationary state is 

!>[/ d 3 rS4>*(r){t + v(r) - ei }^(r)n j¥i (j|j) + cc] 

i 

+ J2n i {i\i+v(v)-e i \i)SU j7ti (j\j). (4) 

i 

This vanishes for a stationary state, but cannot do so 
unless all Lagrange multipliers are included, precluding 
an equivalent Thomas-Fermi theory. 

The argument presented by LS fails to consider the 
full variational expression including (i\v(r) — Ci\i) for each 
occupied orbital. LS Eq.(33) should be written as 



ST S = (<tyi|ei-v(r)|&)(2|2) 

+ (5<h\e2-v(T)\<h)(l\l) + cc, 



(5) 



plus terms involving 5(i\i) that drop out of Eqs.(4). A 
correct derivation must take into account that the nor- 
malization product modifies all terms of Eq.(3), and mod- 
ifies the usual DFT definition p = J2i n iPi- The func- 
tional derivative must be consistent with the Hartree- 
Fock equations implied for Fcrmi-Dirac electrons. Fail- 
ure to include a separate Lagrange multiplier for normal- 
ization of each occupied orbital is inconsistent with the 
constraints needed to specify a Slater determinant as the 
variational trial function. 

The chain rule connecting orbital and density func- 
tional derivatives of the functional T s [p] is a consistency 
condition implied by the Kohn-Sham construction of T s 
using the Schrodinger orbital functional T [7]. In their 
discussion of the chain rule appropriate to DFT, LS in- 
voke an expression, LS Eq.(36), that is undefined unless 
a Frechet functional derivative exists. The integral in 
LS Eq.(36) is inappropriate, since p is an explicit func- 
tion of <f>*. LS Eq.(38) is inconsistent with the orbital 



Euler-Lagrange equations. Such logic can be avoided be- 
cause the ground-state functional, extended to include 
density variations unrestricted by normalization, defines 
Gateaux derivatives [2] such that [7] 
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- =t(pi = {Ei - V}(j)i 



(6) 



A Frechet derivative is implied only if all Ci are equal. 
Alternatively, since these Gateaux derivatives are oper- 
ationally equivalent to the Schrodinger operator t, one 
can define a generalized functional derivative as such a 
linear operator. Using this generalized definition, the 
Kohn-Sham equations can be derived entirely in terms 
of density functional derivatives [8] . 

LS misstate the mathematical implications of the 
fact that the noninteracting Kohn-Sham or Schrodinger 
equations imply by direct construction the existence of 
Gateaux functional derivatives of the extended functional 
T s [7] . Unless all ej are equal, this implies that a Frechet 
functional derivative does not exist. If a Frechet deriva- 
tive is simply assumed, in order to define an extended 
functional of total density, the theory is inconsistent with 
Fermi-Dirac statistics. As a consequence of the exclu- 
sion principle, there is no physically correct definition of 
such an extended functional. However, because Gateaux 
derivatives exist, a unique extended functional of the sub- 
shell densities pi is defined. This suffices to determine the 
Euler-Lagrange equations of the theory [9] . For practical 
applications of DFT, which must assume a parametrized 
functional E xc [p], the implied universal functional F s has 
the same properties as T s . An extended functional F s [p] 
and the Frechet derivative required by TFT do not exist 
in general for more than two electrons, but a unique ex- 
tended functional of the partial densities pi does exist, as 
do the Gateaux derivatives required to determine Euler- 
Lagrange equations. This generalizes TFT to a theory 
consistent with Fermi-Dirac statistics and with electronic 
shell structure [9]. 
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